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1 Introduction 

The study of characteristics of objects of different 
nature carried out in many fields of science and 
engineering is based on measurements of intensity 
of magnetic field created by normal currents, 
flowing in these objects, supercurrents (in 
superconductors) or bound currents of 
magnetization. The results of these measurements 
are the input data for the solution of an inverse 
problem and evaluation of the density of the currents 
which have created the field being considered. 
Among the known methods of numerical research of 
the diffluent currents process in two-dimensional 
area three methods have found practical application: 
finite-element, spatial filterings based on direct and 
a inverse Fourier transform, and probability . 
Without discussing these methods in detail, we shall 
emphasize their essential quality. In each of them 
the unknown is the distribution of the current 
density vectors j. After an evaluation with a selected 
method of their components they are represented by 
arrows, which in size are proportional to the module 
of j. At that the sources with small intensity are lost. 
Such representation of a current is not visually 
enough and in many applications, for example in 
biomagnetic researches, it causes difficulty in 
evaluation of the true location of the current lines. 
Besides, experience of computations shows, the 
precision of the solution of an inverse problem with 
such an approach (analysis of density vectors 
distribution) is essentially reduced as the distance 
between area of measuring of a magnetic field and a 
surface (in many cases flat), containing currents 
increases. 

The author [1] originates other, in comparison with 
traditional, physicomathematical model of sources 
of a magnetic field as a double layer of magnetic 
charges, where the surface density is identical with a 
current function. 

Transition to model of secondary sources in form of 
double layer magnetic charges and appropriate 
integral equation permitted to make a general 
numerical method of reconstruction of surface 
current. This model is suitable for the solution of 
different classes electrophysics problems, for 
calculation and visualization as current lines of a 


two-dimensional currents diffluence. It is 
irrespective of their physical nature (normal, 
supercurrents and currents of a magnetization) and 
from a configuration of the carrier. As the 
elementary source of a field at such formulation 
represents a magnetic moment, the designed 
approach is called Magnetic Moments Method 
(MMM) [1-4]. 

The main difference between MMM and the well- 
known methods consist of the visualization of 
biocurrents as continuous lines. Moreover, a map of 
current density vectors can be constructed as in 
comprehensive methods. 

According to the accepted physical and mathematical 
model of sources of magnetic field the integral 
equation (IE) is constituted and solved, the right 
member of which represents the measured magnetic 
field and the solution is the distribution of its sources 
in investigated plane. 

The level lines (the lines of a constant value) 
obtained after solving this equation represent the 
current lines. 

On a map of current lines the sources of both large 
and small intensity are precisely mapped. The 
current vortexes near on size to calculated area sites 
are mapped irrespective of their amount and 
intensity. Thus the MMM method makes it possible 
to see not only common pattern of biocurrents, but 
also to distinguish its details. 

2 Method 

2.1 Use of method MMM in problems of 

electrophysics. Direct and inverse 
problems 

We shall represent a layer of currents distributed in a 
two-dimensional space (on a surface) S by a double 
layer of magnetic charges by a density x and we 
shall consider originating thus direct and inverse 
problems. Following tradition, we shall formulate 
the direct problem as a computation in the area of 
observation S', induction B or magnetic intensity H, 
as functions from distribution T. 

Superconductors. In the study of the magnetic field 
of superconductors the direct problem should also 
include preliminary computation T on a given 





function H„° , where H n ° is normal to S a component 
of a field strength H° of known (primary) sources. 
The inverse problem is in searching for distribution 
T, corresponding to the measured (secondary) field 

H. 

The model of a simple layer of currents as a double 
layer of magnetic charges can be used for 
visualization not only at the solution of a direct 
problem for flat areas, but also both for coplanar 
superconducting constructions, and for surfaces of a 
general view. Besides it allows to solve a direct 
problem for "massive" superconductors, which 
surface is closed and has an arbitrary shape. This 
relates to calculation of the high frequency currents 
devices as well. 

For the analysis of a magnetic field of 
superconductors the direct and inverse problems are 
formulated as IE of the first kind. The problem of 
visualization of supercurrents, generating this field, 
as current lines is solved. 

Diamagnetic and paramagnetic materials. One way 
of studying weak-magnetic materials and objects is 
that they are placed in homogeneous vertical or 
horizontal a magnetic field. The secondary 
(nonhomogeneous) field generated as a result and 
determined by the magnetization of the material, is 
registered by SQUID detectors in the selected 
(usually flat) area of space. If the cross-section of 
the investigated sample is much less than its other 
two dimensions, with some error, it is possible to 
consider currents of magnetization contained in its 
median plane. Then, having solved the inverse 
problem formulated as the corresponding IE, they 
can be visualized by current lines. 

Biomagnetism. The model of a flat currents 
diffluence is widely used in biomagnetic research. 
The numerical reconstruction of currents from the 
measured magnetic field based on MMM algorithm 
and their depiction by current lines (instead of the 
current density vectors) can improve the resolution 
capability and obviousness of their visualization in 
different organs of the person (heart, brain, liver and 
others). Evidently, that the complete solution of the 
magnetic mapping problem of organs is achieved 
only with 3D problem statement. 

The simplified variant of a inverse problem of a 
biomagnetism is the search for the sources of a 
magnetic field in a plane of parallel flat area of 
measuring. 


2.2 Magnetocardiographic analysis 

One of the methods of medical diagnostics of 
cardiac activity is the analysis of a pattern of 
magnetic field created by biological currents of the 
heart. 

The easiest method of evaluation of cardiac 
performance is the analysis of features of the 
measured magnetic field distribution, in particular 
its vertical (normal in relation to the measurement 
plane) components. Based on this, the attempts are 
undertaken to differentiate separate pathologies. A 
far more informative approach to understanding 
cardiac performance is one grounded in the research 
of distribution of biological currents of the heart 
which, generally speaking, flow throughout its 
volume. 

Due to complexity of this problem stipulated, first of 
all, by its multivaluedness, it is accepted to solve it 
in a simplified form. It is accepted that the sources 
of the field are located in the plane (Fig.l), separated 
from the measurement plane by some distance 
conventionally called the depth of sources ( h). 
Therefore the given inverse problem is formulated 
so: given the measured magnetic field, to find the 
distribution of biological heart currents in a sources 
plane, which with an accepted precision reconstructs 
the measured field. 

In [4] for magnetocardiographic (MCG) analysis the 
MMM method is employed. For its basis the model 
of currents by way of distribution on a plane S of 
magnetic dipoles or current function T was 
accepted. Function ( t) is defined by formula 

j = [Vx^ , n] , (1) 

where j - surface current density , n - unit normal 
to S and Vt^ - surface gradient in flowing point. 
To restore currents in plane under consideration, the 
integral equation of the 1-st kind is to be solved. The 
IE rather T looks like 

\l{M)K{Q,M)dS M =f(Q), (2) 

s 

where (Q, M) designate the point of observation 
and flowing point correspondingly; the kernel K(Q, 
M) corresponds to an output signal gradiometer 
(zero, first or second order) due to a dot magnetic 
dipole disposed in S and having a unit density; / (Q) 
- the right-hand side of the equation is equal to total 
output signals gradiometer, obtained at measurement 
process in nodes of S'. 




Figure 1: Digitization of the plane of measurement 
S' and of the plane S which contains field sources 

The graph of isolines of the calculated function t 
represents a pattern of current lines lying in S. 

For numerical research of the flowing in a heart 
biocurrents simulated by the flat currents the 
program CURLINES realizing MMM- method was 
developed by author [1] . 

As an example the traditional current density map is 
presented in Fig.2 . 
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Figure 2: Image of the heart currents in a source 
plane (S) as the traditional current density map 

The corresponding current lines map is in Fig. 3. 
Both images reflect the distribution of biomagnetic 
sources of heart in a source plane S and are obtained 
by MMM method while running the program 
CURLINES. The size S' and S (Fig.l) in the given 
variant of the program constitute 200 mm by 
200 mm. 


In the map of current density (see Fig.2) each arrow 
displays direction and normalized magnitude of a 
current density in an appropriate point of plane S. 
Current lines (Fig.3) represent lines in each point of 
which the vector of a current density is directed on a 
tangent to this line. A positive value of a current 
function is mapped by light color, negative value - 
by dark color. 



Figure 3: Image of the heart currents in a source 
plane (S) as the current lines map 

The numerous results obtained with the help of the 
MMM method and the program CURLINES permit 
us to generalize that the map of current lines gives 
a more complete and detailed representation of 
currents flowing in the heart in comparison to a map 
of density vectors. 

The visualization of heart currents as current lines 
permits to observe their complete topology, being 
not limited by quantity of concentrated sources 
(current vortex) and by their sizes and intensity 
represent current lines. Moreover using the formula 
(2), it is not difficult to calculate distribution of 
current density also. 

2.3 Magnetographic analysis of brain 

MMM-method can be used for visualization of 
biological currents of a brain with a sufficient 
resolution without restriction of to the flat shape of 
area containing them i.e. other regular surface. 
However there are no key limitations to application 
of the method for visualization of volumetric 
magnetic field sources [3], 








3 Results 

The study of the resolution capability of the 
magnetic moments method (MMM) was conducted 
with localization of concentrated current contours 
commensurable to a grid step. It was discovered that 
with increase of the degree of digitization the 
resolution capability at visualization of sources of 
small size increases. The location of the 
reconstructed sources corresponds to their given 
location in the solution of the direct problem. 

The numerical experiments on visualization of 
currents of the heart in a calculated plane were 
conducted. Thereby the right member of IE was 
determined not by the calculated magnetic field 
found by the solution of a direct problem, but by 
measurement data in clinical conditions. 

We have to emphasize that biological currents of the 
heart actually have volumetric character and their 
image in a calculated plane is approximated to a 
very large degree. Despite of it, the obtained planar 
current distribution is such that the magnetic field 
calculated on it coincides close enough with the 
measured field. The map of current lines allows to 
analyze the course of a cardiocycle in all of its 
phases, which could not be done by traditional 
method. 

Besides the map of currents on the found 
distribution of a density of double layer magnetic 
charges it is easy to construct the usual map of 
current density vectors, and in this case a high 
enough accuracy of evaluations is preserved for 
vectors of small magnitudes as well. 

The method allows to distinguish close enough 
posed moments of magnet (curls of a current) with is 
opposite directional vertical axis (see [4]). From it 
by the help it is possible to renovate horizontally 
directional moments of magnet (curls) in a 
horizontal plane of sources (see [3]). 

Further work on perfecting of the method should 
consist of increase of the resolution capability and 
the degree of adequacy of flat current map to actual 
distribution of volumetric currents of the heart. 

4 Discussion 

The resolution of a method MMM depends 
somewhat on a relation of intensities of sources of a 
magnetic field, and also from their mutual location 
in space. 

With increase of distance of a source from a plane of 
measuring it the magnetic field is attenuated. 
Together with it blur a fashion of a source in a 
plane of sources. 


To receive an acceptable resolution of a currents 
pattern, is necessary to introduce enough high 
degree of task digitazaiton. It leads to its bad 
conditionality, which will increase with a raise of 
an appropriate matrix order. Restoring of current 
distributions with desired accuracy need further 
perfecting of the numerical methods for solution of 
formulated above task. 
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